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ABSTRACT 

Galaxies on the largest scales of the Universe are observed to be embedded 
in the filamentary cosmic web which is shaped by the nonlinear tidal field. As 
an efficient tool to quantitatively describe the statistics of this cosmic web, we 
present the anisotropic two-point correlation functions of the nonlinear traceless 
tidal field in the principal-axis frame, which are measured using numerical data 
from an N-body simulation. We show that both of the nonlinear density and 
traceless tidal fields are more strongly correlated along the directions perpen- 
dicular to the eigenvectors associated with the largest eigenvalues of the local 
tidal field. The correlation length scale of the traceless tidal field is found to be 
~ 20 ft, -1 Mpc, which is much larger than that of the density field ~ 5 h~ 1 Mpc. 
We also provide analytic fitting formulae for the anisotropic correlation functions 
of the traceless tidal field, which turn out to be in excellent agreement with the 
numerical results. We expect that our numerical results and analytic formula are 
useful to disentangle cosmological information from the filamentary network of 
the large-scale structures. 

Subject headings: cosmology: theory — large-scale structure of universe 

1. INTRODUCTION 

On the largest scales, the spatial distribution of galaxies exhibits the pattern of a fil- 
amentary network which is commonly referred to as the "cosmic web". N-body simula- 
tions have long predicted the web-like pattern in the clustering of matter and galaxies (e.g., 
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Shandarin fe ZelT3ovichlll986[) . recent large galaxy surveys have confirmed its presence (e.g., 
Colless et al.ll200ll ). and theoretical endeavors have been able to expla in the naturalness of 



its occurrence in the context of the cold dark matter (CDM) paradigm (IBond fc Myerslll996 



Bond et al.lll996l ). It is now well accepted that the cosmic web is a real, natural, and common 
phenomenon. 

Nonetheless, the most crucial part has been missing for the comprehension of the cosmic 
web - an optimal statistical tool to quantify it. The standard N-point statistics of course 
can be used in principle to quantify it, but they cannot capture the essence of the cosmic 
web that it is a nonlinear manifestation of the primordial tidal interactions. It is because 
in N-point statistics what is emphasized is the role of gravity in establishing the high-order 
correlations of the density field. The cosmic web is a nonlinear manifestation of the spatial 
coherence of the tidal field (IBond et al.lll996l ). Thus, to quantify the cosmic web and capture 
its essence, it is necessary to have a statistical tool by which the role and effect of the tidal 
field and its spatial correlation are highlighted. 

Patchy efforts have so far been made to account for the effect of the tid al field on struc- 
ture formation. For instance, the failure of the Press-Schecther mass function ( Press fc Schechter 



1974 ) has been attributed to the deviation of the gravita t ional process from spherical dynam- 



ics caused by the effect o f the tidal forces (lMonacolll995l ; lAudit et al.lll997l ; iLee fc Shandarin 



19981 ; ISheth et al.l 1200 ll ). The galaxy angular momentum has been ascribed to the ini- 



tial tidal interactions with the surrounding matter distribution (|Doroshkevichl 119701 ; iWhite 
1984J : ICatelan fe TheunJl996l : ILee fc Penlkoool : IPorciani et "alll2002al lbh. The intrinsic galaxy 
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2001; 
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2005; 


Lee &; Pen 



20071 . 120081 ). The alignment and ellipticity correlation functions have been developed to quan- 



tify the preferential orien tations of galaxies relative to the surrounding large-scale structure 



( jFaltenbacher et al.ll2009l ). The cosmic web, however, requires not a patchy but a unified 



way in which the multiple aspects of the nonlinear tidal field can be quantified. 

The difficulty in quantifying the effect of the tidal field lies in th e fact that the cor rela- 
tions of the tidal field are h ighly anisotropic even in the linear regime. ILee &; Pen! (120011 ) and 
Catelan fc Porcianil (120011 ) have for the first time derived analytically the anisotropic two- 
point correlations of the linear tidal field. But, their results were all expressed in the fixed 
Cartesian coordinate frame. To fully appreciate the anisotropic nature of the correlations of 
the tidal field, it is desirable to find the correlations of the tidal field in the principal axis 
frame, since the filamentary cosmic web i s generated by th e anisotropic collapse of matter 
along the principal axes of the tidal field (IBond et al.lll996l ). Unfortunately, it is extremely 
difficult to derive analytically the correlations of even the linear tidal field in the principal 
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axis frame which rotates from point to point. 

The goal of this Paper is to measure numerically the anisotropic two-point correlations 
of the nonlinear tidal field in the principal-axis frame and to provide an analytic fitting 
formula as an efficient statistical tool for the description of the cosmic web. This paper is 
organized as follows. In §2, we provide a brief description of the numerical data. In §3, we 
present the numerical results and analytic fitting formula as well. In §4, we summarize the 
results and discuss the implications. 



2. DATA 



We use the N-body results derived by lHahn et al.l (120071 ) using the tree-PM code GADGET 



2 (ISpringel et al.l 120051 ) . The simulations followed the evolution of 512 3 particles starting 
from redshift z = 52.4 to z = in a periodic box of linear size 180/i -1 Mpc. A flat ACDM 
cosmology is assumed with key parameters given as fl m = 0.25 (matter density), = 0.75 
(vacuum energy density ), fl& = 0.045 (baryon density), cr 8 = 0.9 (power spectrum ampli- 
tude), Hq = 0.73 (Hubble constant) and n s — 1 (spectral index). 



Hahn et al.l (120071 ) constructed the overdensity field 5 (dimensionless density contrast) 
on 1024 3 grids with the help of the Cloud-In-Cell interpolation method and smoothed it by 
a Gaussian kernel with a filtering radius R s . The peculiar gravitational potential field 
was also derived from the overdensity field by solving the Poisson's equation V 2 oc 5. In 
simulations, the overdensity field is converted to the Fourier space with the help of a fast 
Fourier transform (FFT). The potential field is obtained as oc k~ 2 5 where k represents 
the magnitude of the wave vector in the Fourier space. For a detailed descri ption of the 



N-bod y simulations and the derivation of the density and potential fields, see lHahn et al. 



(120071 ). For the smoothing of the density (and potential) fields, we use the small smoothing 
scale R s = 0.5/i _1 Mpc to minimize the effect of the smoothing on the determination of the 
correlations of the tidal fields. 

The tidal tensor (Tij) is defined as the Hessian of the peculiar gravitational potential 
Tij = didj <p . In Fourier space, it is written as = kikj . Using the FFT again, we 
calculate the tidal tensor at each of the randomly selected 250000 grid points. The sum of 
the diagonal components of (T^) (i.e., the trace of the tidal field) equals the dimensionless 
density contrast, 5. To sort out the true tidal effect from the gravitational effect, we cal- 
culate the traceless tidal field (Tij) as = — (5/3) where (iy) denotes the identity 
tensor. Diagonalizing at each of the 250000 sampling points, we find its three eigenvalues, 
{Ai, A 2 , A 3 } (in a decreasing order, i.e., Ai > A 2 > A 3 ), and the corresponding normalized 
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eigenvectors with unit magnitude, {ei, e 2 , 63}. Regardless of whether a given volume ele- 
ment is expanding or contracting (negative or positive S) the traceless tidal tensor gives 
the deformation of that volume element (relative to spherical expansion or contraction). By 
construction, we have Ai > and A 3 < 0. The second largest eigenvalue A 2 can take on 
either sign. The maximal (relative) compression of a given volume element by gravity occurs 
along the direction parallel to ej associated with Ai. Figure Q] plots the probability density 
distributions of Ai, A 2 and A 3 as solid, dashed and dotted lines, respectively. The probability 
density distribution of Ai and A 3 is positively and negatively skewed, respectively while that 
of A2 is almost symmetric about zero. Note that the two distributions, p(Ai) and p(A 3 ), are 
also almost symmetric to each other around zero: p(Ai) ~ p(— A3). The distribution p(A 2 ) is 
much narrower than p(Ai). 

In the following section, we measure the two-point correlation functions of Aj, A 2 , A 3 
and S in the principal-axis frame of the traceless tidal field and study their characteristic 
behaviors. 

3. NUMERICAL ANALYSIS AND FITTING FORMULA 

For each pair of points located at x and x + r, we express the separation vector r in the 
frame of the principal axes of Tij(x) spanned by {e 1; e 2 , e 3 } with orienting the z-axis in 
the direction of ei (the first principal axis of the local tidal field). Let (r, 6, 0) represent the 
spherical polar coordinate of r in the principal-axis frame. For every grid pair, we calculate 
the product, Aj(x)Aj(x + r) (for i = 1, 2, 3). Then, we take the average of the product over 
the points whose separation vectors belong to the differential bin, [r, r + dr] , to determine 
the correlation of A, function of (r, 6, <fi) as 

^(r,^,0) = (A 4 (x)-A i (r + x)). (1) 

Figure [2] shows the isotropic two-point correlations ^\i{r) of the three eigenvalues of 
the traceless tidal field as a function of r (solid, dashed and dotted line for Ai, A 2 and A3, 
respectively) at z = 0. This isotropic correlation function £,\i(r) is obtained by taking the 
average Ai(x)Aj(x + r) over the points whose separation distances belong to the differential 
distance bin, [r,r + dr]. The tendency of £,\i(r) being flat for r < 2 /i -1 Mpc is due to the effect 
of the smoothing of the tidal field. The isotropic two-point correlation of S is also plotted for 
comparison (thin solid line). As it can be seen, the correlations of Ai and A 3 decrease with r 
much more slowly than that of 5. In other words, the largest and smallest eigenvalues of the 
traceless tidal field are coherent over larger scales than the density field 5. Let the correlation 
length scale l c correspond to the distance at which the correlation strength becomes an order 



5 



of magnitude lower than its value at r = 0. We find that l c ~ 20/i~ 1 Mpc for Ai and A3, while 
5 and A2 has l c ~ 3h~ 1 Mpc. The smallest eigenvalue A3 has the highest amplitude at all 
distances, while the second largest eigenvalue A 2 has the lowest amplitude at r < lO/i^Mpc. 
The correlation amplitude of 5 is comparable to that of Ai at small distance r < 2/i~ 1 Mpc 
but decreases much more rapidly with r than Ai, reaching the lowest value at r > 10/i _1 Mpc. 

For the rescaled isotropic correlation ^Ai( r )/^Ai(0), we find the following analytic formula: 



where tq and n are two fitting parameters. By comparing the numerical results with the 
above analytic formula, we empirically determine the best-fit values of r and n, which are 
listed in Table 1. Figure [3] shows the rescaled isotropic correlations of Ai (left), A2 (middle) 
and A3 (right) at z = 0. In each panel the dots and the solid line represent the numerical 
results and the analytic model with best-fit parameters, respectively. The error at each bin 
is calculated as one standard deviation in the measurement of the mean value. As it can be 
seen, the analytic models (eqs.[2]-[3]) fit the numerical results very well for r > 2 h^Mpc. 
On smaller distances r < 2/?r 1 Mpc, the numerical results become flatter than the analytic 
curve due to the smoothing effect. 

Figure H] shows the anisotropic two-point correlations of Ai (left), A2 (middle) and A3 
(right) as a function of r and their variation with 9 at z = 0. These anisotropic correlations 
£,\i(r, 9) are obtained by taking the average of Aj(x)Aj(x+r) over the points whose separation 
distances belong to the differential distance bin [r, r + dr] and at which the angles between 
the separation vectors and the first principal axes are in the differential range of [9,9 + d9\. 
In each panel the thick solid, dashed, dotted, longdashed, dot-dashed, and dot-longdashed 
line correspond to the range 0° < 9 < 15°, 15° < 9 < 30°, 30° < 9 < 45°, 45° < 9 < 60°, 
60° < 9 < 75° and 75° < 9 < 90°, respectively. The fluctuations of ^(f, #) in the small-r 
section (r < 2 /t _1 Mpc) are due to the numerical noise. The anisotropic two-point correlation 
of 5 is also plotted for comparison as thin lines in each panel. As it can be seen, both 
correlations of A and 5 increase with 9 in the principal-axis frame at all distances. That is, 
the correlations become stronger toward the directions normal to the first principal axes of 
the tidal field. This behavior is consistent with the scenario that matter compress maximally 
along the first principal axes of the tidal field and thus is more likely to lie in the plane normal 
to the first principal axis of the tidal field. The variation of A2 with 9 is stronger than that of 
Ai and A3. Note also that the degree of variation of ^Ai( r , 0) with 9 decreases with r. In other 
words, the correlations tend to become less anisotropic at large distances (r > 10 /i _1 Mpc). 




(2) 



(3) 
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We also obtain £^(0) by taking the average of A;(x)A;(x+r) over the points at which the 
angles between the first principal axis and the separation vectors belong to the differential 
angle bin, [9, 9 + dd]. The following analytic formulae are found to fit the numerical results 
well: 

— = 1 -^7, for A X ,A 3 , (4) 



£(0) 1-Asin 4 9 
f(0) ~" 1-Asm 2 9'' 



for A 2 (5) 



where A is an adjustable parameter whose best-fit value is determined via the \ 2 statistics 
(see Table 1). Figure 0] shows the anisotropic correlation function of Ai (right), A 2 (middle) 
and A3 (left) as a function of 9 in the principal axis frame of the local tidal field. In each 
panel the dots represent the numerical results while the solid line is the analytic fitting 
formula. As it can be seen, the analytic models fit the numerical results very well for each 
case. To improve readibility we do not show errorbars which are negligibly small anyway. 
These results show clearly how anisotropic the correlations of the eigenvalues of the traceless 
tidal field are with respect to the first principal axis. The difference between the correlations 
at 9 = and 9 = 90° reaches approximately 60% for Ai and A3, while it increases up to 70% 
for X 2 . 



4. SUMMARY AND DISCUSSION 

The web-like pattern in the large-scale matter distribution is an imprint produced by the 
effect of the nonlinear tidal field in the late universe. While the trace part of the tidal field is 
responsible for the collapse of matter, its traceless part deforms the cosmic volume elements, 
resisting the overall spherical expansion or contraction. In consequence the traceless tidal 
field induces anisotropy in the large-scale correlations of matter distribution. An efficient 
statistical tool has long been looked for to describe quantitatively this cosmic web and to 
retrieve cosmological information encoded in the cosmic web. 

We have determined numerically the two-point correlations of the three eigenvalues 
of the nonlinear traceless tidal field defined as the Hessian of the gravitational potential 
in the frame of the principal axes of the tidal field. The numerical findings indicate that 
the correlation functions of the traceless tidal field and the density field are all anisotropic 
relative to the principal axes. Their correlations having much larger correlation length scales 
than that of the density field, increase along the directions normal to the first principal 
axes of the tidal field. The analytic fitting formula for the correlation functions of the three 
eigenvalues of the tidal field are determined and found to fit the numerical results well. 
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Our numerical and analytical results are thus able to provide a physical description of the 
large-scale filamentary structure. 

It will be interesting to know how the anisotropic two-point correlations of thetraceless 
tidal field depend on the background cosmology. Our future work is in this direction. 

J.L. acknowledges the financial support from the Korea Science and Engineering Founda- 
tion (KOSEF) grant funded by the Korean Government (MOST, NO. RO 1-2007-000-10246- 
0). O.H. acknowledges support from the Swiss National Science Foundation. All simulations 
were performed on the Gonzales cluster at ETH Zurich, Switzerland. 
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Fig. 1. — Probability density distributions of the three eigenvalues, Ai, A 2 and A3 (solid, 
dotted and dashed line, respectively) of the traceless tidal field smoothed on the scale of 
0.5/i _1 Mpc. 



- 11 - 



T 




1 



1 10 

r [h _1 Mpc] 

Fig. 2. — Two-point isotropic correlation functions of the three eigenvalues of the traceless 
tidal field, Ai, A 2 and A 3 (solid, dashed and dotted line, respectively) smoothed on the scale 
of 0.5/i -1 Mpc. The two-point isotropic correlations of 5 is also plotted (thin solid line) for 
comparison. 
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Fig. 3. — Rescaled isotropic correlation functions of Ai (right), A 2 (middle) and A 3 (left) in 
the principal-axis frame of the local tidal field. In each panel the dots represent the numerical 
results while the solid line is the analytic fitting formula. The errors represent one standard 
deviation in the measurement of the mean values. 
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Fig. 4. — Anisotropic correlation functions of Ai (right), A2 (middle) and A3 (left) as a 
function of r and their variations with polar angle 9 in the principal-axis frame of the local 
tidal field (thick solid lines). The polar angle represents the angle between the position 
vector and the first eigenvector of the local tidal field. In each panel the solid, dotted, dashed, 
long-dashed, dot-dashed, and dot-longdahsed line correspond to the range, 0° < 9 < 15°, 
15° < 9 < 30°, 30° < 9 < 45°, 45° < 9 < 60°, 60° < 9 < 75° and 75° < 9 < 90°, respectively. 
The anisotropic correlation functions of 5 are also shown for comparison (thin solid lines). 
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Fig. 5. — Anisotropic correlation function of Ai (right), A2 (middle) and A3 (left) as a function 
of 9 in the principal axis frame of the local tidal field. In each panel the dots represent the 
numerical results while the solid line is the analytic fitting formula. 
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Table 1. Best-fit values of the parameters in the fitting formula for the two-point 
anisotropic correlation functions of the three eigenvalues of the nonlinear traceless tidal 

field. 



parameters 


Ai 


A 2 


A 3 


n 


1.28 


2.76 


1.65 


ro 


3.03 


12.58 


21.51 


A 


0.37 


0.41 


0.38 



